Brownian forces in sheared granular matter 
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We present results from a series of experiments on a granular medium sheared in a Couette 
geometry and show that their statistical properties can be computed in a quantitative way from 
the assumption that the resultant from the set of forces acting in the system performs a Brownian 
motion. The same assumption has been utilised, with success, to describe other phenomena, such as 
the Barkhausen effect in ferromagnets, and so the scheme suggests itself as a more general description 
of a wider class of driven instabilities. 
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The shear response of granular media has been in- 
vestigated for more than a century by different scien- 
tific communities jl| , ranging from soil mechanics |2| and 
earth sciences 0, Q to physics 0, 0|. This fact is re- 
flected by the variety of experimental settings and the- 
oretical methodologies applied to the problem, such as 
rate and state friction laws 0, or direct particle sim- 
ulations Under a slow loading rate, the shear re- 
sponse of gran ular media typically displays large fluctu- 
ations [lfj, Ull , due to the rearrangements of the inter- 
nal force chains [T3 . fl3| along which stress propagates in 
granular materials. 

Motivated by the practical concern of obtaining a me- 
chanical description to be used in applications, most of 
the past theoretical activity has been devoted to the anal- 
ysis of averaged properties, in the search for constitutive 
macroscopic laws. A detailed analysis of slip statistics 
is, however, of primary importance if one wishes to un- 
derstand the mechanisms which are at the base of the 
wide fluctuations observed. They are the distinguish- 
ing feature of " crackling noise" , which denotes an inter- 
mittent activity with widely fluctuating amplitudes |l4j . 
Crackling noise is observed in different contexts such as 
magnetic materials 0, Htl. f erroelectrics. type II super- 
conductors [13, fracturefta. Il9| and plasticity pel l2l| 
to name just a few. In addition it shares analogies with 
earthquake phenomenology, where the statistics typically 
displays wide fluctuations and recorded sizes of seismic 
events span several decades [22|, 1.23] . Common aspects 
observed in such apparently different systems have sug- 
gested the existence of a deeper correspondence in the 
underlying physical properties 0], and one may hope 
that any general insight or understanding gained in these 
phenomena may guide research in other fields. Among 
the schemes proposed in this context, perhaps the most 
successful are the scaling-renormalization group |14| and 
Self-Organized Criticality |24|. Within these frameworks 
many phenomena have been qualitatively explained in 
terms of a set of logical, coherent and simple rules, but it 
is generally difficult to adapt them to specific situations 
and obtain quantitative predictions. In this letter, start- 
ing from a laboratory experiment on sheared granular 



FIG. 1: The signals measured in experiments: a) instanta- 
neous velocity of the disk; b) instantaneous torque exerted by 
the spring; c) definition of duration and size of a slip event. 



matter we propose a model which quantitatively repro- 
duces the observations and, at the same time, is based 
on a simple but intuitive hypothesis which sheds light on 
the similarities between different phenomena. 

Our experimental set-up consists of a Couette cell 
made of a circular channel containing mono-disperse 
2mm glass beads. An annular plate is driven over the 
top surface of the channel by a motor through a torsion 
spring with stiffness k. In most of the experiments we em- 
ploy a spring of stiffness k = 0.36N, but we also consider 
different values (i.e. k = 0.12N and 0.9N). To ensure a 
granular shearing plane, the annular top plate has a layer 
of beads glued to its lower surface, though, to avoid indi- 
vidual grains jamming the system at the boundary, this 
layer of beads does not extend the full width of the chan- 
nel. The system is initialized before each experiment by 
pouring the medium into the channel, then the system is 
run at a slow velocity V (~ 0.01 rad/s) for a long time, 
in order to approach a stationary state. We have ob- 
served that the mean torque approaches a steady-state 
after the order of 100 revolutions. Our measurements 
from the device consist of the angular position of the 
top plate and the deflection of the torsion spring, sam- 
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FIG. 2: The experimentally measured distribution of disk 
velocities for different values of the driving angular velocity 
V. The lines are obtained from the numerical integration of 
Eqs. 11131 using parameters value extracted from the exper- 
iments. The distributions for different drives are vertically 
shifted for clarity. 

pled at 1 kHz with an error ~ 10~ 4 rad. From this, 
with knowledge of the torsion spring constant and the 
moments of inertia of the plate, we obtain angular ve- 
locity and acceleration, instantaneous torque and work 
done by the spring. Varying the driving angular velocity 
V in the range 10~ 4 -10 _1 rad/s, the system displays a 
the stick-slip regime: the disk is still as the spring winds 
up and the torque on the plate increases until the plate 
slips. The equation of motion of the disk reads: 

19 = k(Vt -6)-F, (1) 

where I is the moment of inertia of the system, k is the 
stiffness of the torsional spring, and V is the angular ve- 
locity of the driving motor. The first term on the right- 
hand side is the force exerted by the spring. The motion 
of the plate is counteracted by the " friction" F exerted 
by the medium. Is this last term that encompass the 
whole complexity of the physical stick-slip phenomena. 

The typical instantaneous velocity signal displays a 
very erratic behavior with pulses of widely fluctuating 
magnitudes (see Fig-QJi), and its probability density dis- 
tribution spans more than two decades. The presence of 
power law statistics in the stick-slip phase of a sheared 
granular was first reported some years ago in slightly 
different experiments [Tlj The present experiments show 
that the general shape of the distributions depends on the 
shear rate, at odds with the cutoff position (see Fig. |2). 
One could attempt to fit the data in Fig. [3 with a power 
law to obtain an exponent, but the limited scaling regime 
would not yield a reliable result. In order to get a more 
complete characterization of the fluctuations and their 
temporal scales, we analyze duration T and a size S as- 
sociated to each pulse, as is commonly done in crackling 
noise experiments in this case S is equal to the an- 
gular slip of the plate during the pulse (see Fig.^). The 
distributions of T and S, measured for different values 
of the driving angular velocity V, are reported in Fig. [3] 
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FIG. 3: Slip sizes probability distribution, for different values 
of the driving angular velocity. The lines are the results from 
the numerical integration of Eqs. 11 1 lot . The distributions for 
different drives are vertically shifted for clarity. 

and Fig. 0] They display a complex shape, in which an 
initial decay is followed by a characteristic peak at larger 
scales, indicading broad fluctuations over many scales. 
Note that the position of the peaks does not depend on 
the driving velocity V. 

This behavior is entirely contained within F of Eq. QJ, 
the instantaneous frictional torque. Because of the disor- 
dered arrangement of grains, forces propagate within the 
medium via a complicated network in which both elas- 
tic and frictional interaction play fundamental roles |25j . 
Moreover the network is constantly evolving in a pseudo- 
plastic way [2(j , and so a detailed force description rep- 
resents a formidable problem with a very large number 
of interdependent degrees of freedom. The friction term 
F is, in principle, a random function which may depend 
on a variety of state variables The minimal formu- 
lation retains only the dependence on the position and 
velocity of the disk, i.e. we assume that F — F(9,9). 
Note that the dependence on the time arises only from 
9 and 9, since when the disk is stuck the force does not 
change. As an additional assumption, we split F into 
two additive terms: In the first of them, the dependence 
on velocity is confined to a deterministic law, analogous 
to the viscous reaction of a fluid, while the second term 
is independent of the velocity and introduces the fluctu- 
ating part: F(9, 9) = F(9) + F f (9). 

In Fig.[SJi we report the value of the friction torque as 
a function of the instantaneous angular velocity of the 
plate. Once fluctuations have been averaged out, a clear 
velocity dependence is seen, with strain-rate weakening 
at low velocity followed by strengthening at high veloc- 
ity. These features, not present in the viscous reaction 
of a Newtonian fluid, are responsible for the stick-slip 
instability. An adequate fit of the experimental data is 
provided by a rate dependent friction law of the type 

F(9) = F Q + 1 (9-2v Q ln(l + 9/v Q )), (2) 

where Fq is the average static friction torque, vq cor- 
responds to the minimum in the average torque and 7 is 
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FIG. 4: Slip durations probability distribution, for different 
values of the driving angular velocity. The lines are the re- 
sults from the numerical integration of Eqs. 11 lot . The dis- 
tributions for different drives are vertically shifted for clarity. 
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FIG. 5: Frictional torque: a) average velocity dependence; b) 
power spectrum of the random frictional torque. 



the high velocity damping (see Fig. [5]}) . 

It is interesting to note that an equation similar to 
Eq. (J2J) is used to describe the friction law observed in 
dry solid-on-solid friction |27| (steady sliding regime). 
Previous work have shown that the friction force from 
granular materials is not a single valued function of the 
velocity during a single slip event |||. Our analysis does 
not contradict this observation, since the result displayed 
in Fig.[SJi describes the average over different slip events. 

The fluctuating part of the friction torque is, as dis- 
cussed above, a direct consequence of the disordered 
structure of the force chain network present in the gran- 
ular medium. As the disk slips by a small angle SO, one 
expects that the friction torque is increased or decreased 
by a small random amount SFf , reflecting a limited rear- 
rangement in the granular structure. On the other hand, 
under subsequent or large slips, rearrangement of grains 
will be more complete, and the random torque will even- 
tually decorrelate. The simplest mechanism to account 
for these qualitative features is a confined Brownian pro- 
cess: 



dFf_ 
dO 



(3) 



where 77 is an uncorrelated noise term extracted from 
a Gaussian distribution, i.e. (?7(#)^(#')} = D5(6 — 6'), 
where D is the noise variance. It is important to note 
that the independent variable here is the angular position 
9 instead of the time t, as always the case for quenched 
athermal disordered systems. 

The parameter a determines the inverse correlation 
length and can be estimated by the power spectrum S(w) 
of the process. In fact this latter is related to the cor- 
relation function by the Wiener-Kintchine theorem. For 
the process above the power spectrum is given by 



S(k) 



d6F f {0) exp(-i(9fc) 
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(4) 



The power spectrum measured in the experiment is 



shown in Fig. 03- It is seen that despite the simplic- 
ity of Eq.©, the related Lorentzian spectrum of Eq.QJ 
compares quite well with the actual spectrum and allows 
estimation of the inverse correlation length a. Once the 
parameters a, D, Fq, 7 and vq have been extracted from 
the data, the equation of motion JJ| is integrated numer- 
ically, and a stochastic velocity signal is obtained, to be 
compared with the experiments. It should be emphasized 
that the the fitted parameters display quite stable values 
upon repeating the experiments and varying the driving 
rate. Typical values at low V are a « 15 rad^" 1 , D w 0.01 
(N 2 /rad), F » 0.55 Nm, 7 w 0.20 N m s /rad, v Q « 0.03 
rad/s. At high velocities (V > 0.1 rad/s) we observe a 
significant variation of these values, indicating the onset 
of the fluidization transition [28j . 

The model reproduces the complex phenomenology 
of the granular dynamics, despite being based on few 
mean-value parameters whose values are directly mea- 
sured from the system under study. The only unknown 
parameter here is the moment of inertia of the system 
/. In fact, considering I as the bare moment of inertia 
of the disk Iq , one neglects the possibility that a number 
of grains could move together with the plate, increasing 
the actual moment of inertia of the system. On the con- 
trary, this behavior is expected for the granular system, 
due to formation of shear bands j2(| . Accordingly, a best 
quantitative agreement between model and experiments 
is obtained considering I = 1.5/o, where I = 0.02 Kg 
m 2 , corresponding to few granular layers moving together 
with the disk. 

Examples of the distributions resulting from the model 
are compared to experimental data in Fig. [3 Fig. [3] and 
Fig.Q] We notice that the best agreement is found at low 
driving. Possible reasons for the discrepancies observed 
at high drives are related to the presence of the fluidiza- 
tion transition and the intrinsic approximations of our 
model. In particular, the power spectrum is not exactly 
Lorentzian and the model does not account for dilatancy 
effects and complex shear band dynamics. Nevertheless, 
the simplicity of the proposed approach is quite attrac- 
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slip 


e 


domain wall jump 


shear rate 


kV 


magnetic field rate 


spring constant 


k 


demagnetizing factor 


frictional damping 


7 


eddy current damping 


random friction 


Ff 


pinning forces 


moment of inertia 


I 


domain wall mass 



TABLE I: Table of correspondences between physical quan- 
tities for the shear response of a granular medium and the 
domain wall motion in a ferromagnetic material. 



tive and suggests the tantalizing idea that phenomena ob- 
served in different contexts could be effectively described 
by similar general processes. 

This idea is supported by the observation that a sim- 
ilar process emerges in a very different context from the 
present one, specifically that of ferromagnetic materials. 
The equation of motion we adopt to describe granular 
media is similar to the one employed for domain wall 
motion [3(| • The latter reproduces the statistical prop- 
erties of the Barkhausen noise recorded during hysteresis 
in soft magnetic materials. In particular, the granular 
equation reduces to the domain wall model in the limits 
V >> vo and 7 >> I ■ V, as is appropriate for metallic 
ferromagnets where eddy current damping cancels iner- 
tial effects. The random structure of the friction force in 
the case of domain walls is due to the inhomogeneities 
present in the material and the correlation properties are 
believed to arise from a summation of the pinning forces 



along the domain wall [3l| , while here they are due to the 
formation and destruction of force chains. A summary 
of the mapping between granular matter and magnetic 
domain walls is reported in Table [I] 

The domain wall model is exactly solvable and predicts 
a power law dependence for the distribution of velocities, 
slip sizes and durations |30| ; the corresponding exponents 
decrease linearly with the driving rate. The power laws 
are terminated by an exponential cutoff inversely pro- 
portional to the demagnetizing factor [31j . 

Our granular model presents an additional feature: the 
distribution tails display a peak that masks the power law 
decay. This is due to the combined effect of inertia and 
the shear rate weakening frictional torque. This idea is 
confirmed by the fact that the peak in the duration dis- 
tribution scales as T p ~ y/I/k. This behavior is simply 
recovered considering the deterministic motion obtained 
by suppressing friction fluctuations in Eq. QJ, and it has 
been confirmed in our experiments. Future experiments, 
designed to reduce inertial effects and shear rate weaken- 
ing, could allow the observation of the power law decays, 
whose exponents are predicted by our model. 

The results presented here, supported by the results 
for the Barkhausen effect in ferromagnets [3(J], indicate 
that the Brownian scheme (eq. [3J for the description of 
random forces evolution in dissipative systems driven by 
an external force may constitute an effective approach 
of general validity. Similar analysis of different systems 
exhibiting " crackling noise" |14| ] will constitute a proving 
ground for these ideas. 
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